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Abstract. The Crab-like SNR G21.5-0.9 was observed in the X-ray band (0.5-10 keV) by the XMM-Newton satellite for 
over 100 ks. The large effective area of the EPIC instrument has allowed us to to perform a deep search for pulsations from 
the central core of G21.5-0.9 . No pulsations were found with upper limits on the pulsed fraction between 7.5 % and 40 % 
(depending on frequency and energy range). 

Key words. ISM: individual (G21.5-0.9) - supernova remnants - X-rays: ISM 


1. Introduction 

Almost 5% of the ~ 225 known supernova remnants (SNRs) in 
our galaxy are classified as “Crab-like” or “plerionic” (Green 
2000). From the spatial point of view, they are characterized by 
compact, centre-filled radio and X-ray morphology (Weiler & 
Panagia 1978). These SNRs show featureless power-law spec¬ 
tra, with a relatively flat spectral index in the radio regime 
(ar ^0.0 - 0.3) and a steeper one at shorter wavelengths, which 
are typical of synchrotron processes. 

These morphological and spectral characteristics are ex¬ 
plained by the presence of a central pulsar, which injects high 
energy electrons that suffer synchrotron radiation losses as 
they diffuse through the surrounding magnetic field (see, e.g., 
Reynolds & Chanan 1984). 

G21.5-0.9 shows many characteristics of the Crab-like 
remnants. Both in the radio band (Morsi & Reich 1987) and at 
X-ray energies (Becker & Szymkowiak 1981) its emission is 
centrally peaked. Evidence for a non-thermal X-ray spectrum 
was also indicated by GINGA observations (Asaoka & Koyama 
1990). 

The available measurements of the neutral hydrogen ab¬ 
sorption give a distance of ~ 4.8 kpc (Becker & Szymkowiak 
1981). The radio luminosity of G21.5-0.9 is ~ 1.8 x 10^^ d\ erg 
s“^ (Morsi & Reich 1987), i.e. a factor ~ 9 smaller than that 
of the Crab, but its X-ray luminosity is a factor ^100 less, 
therefore the Lx/Lr ratio is significantly lower. 

X-ray observations performed with the Chandra satellite 
detected a compact central core of ^ 2" in size, at the cen¬ 
ter of the more extended synchrotron nebula of ~ 30" radius 
(Slane et al. 2000). The central core, which is spatially re¬ 
solved, most likely marks the position of the pulsar powering 
G21.5-0.9 . Also a fainter, more extended “halo” (radius ~ 2') 
was detected with Chandra. This was tentatively interpreted 
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as the outer “shell” formed by the expanding ejecta and the 
passage of the supernova-driven blastwave (Slane et al. 2000). 
However, a recent XMM-Newton observation shows that also 
the halo has a non-thermal spectrum; it is probably a low sur¬ 
face brightness extension of the plerionic nebula (Warwick et 
al. 2001). However, this interpretation is not supported by radio 
data, since no significant radio emission has yet been detected 
at such a large distance from the source core (Bock & Wright 
2001 ). 

Up to now no pulsed emission has been detected in the ra¬ 
dio or X-ray energy range from the putative neutron star at the 
center of G21.5-0.9 (Frail & Moffet 1993; Kaspi et al. 1996; 
Biggs & Lyne 1996; Slane et al. 2000). In this paper we report 
the results of a sensitive timing analysis on data provided by 
fonr XMM-Newton observations. 


2. Observations and data analysis 

During April 2000, G21.5-0.9 was observed four times by the 
XMM-Newton mission as one of the calibration targets. A log 
of the observations is given in Table 1. The source was on-axis 
in the first observation and ~ 10' off-axis in the following ones; 
in each case the accumulated exposure time was ~ 30 ks. The 
results reported by Warwick et al. (2001) were based only on 
the on-axis observation. 

In all the observations the source was imaged by the three 
focal plane CCD cameras (MOSl, MOS2 and PN) of the EPIC 
instrument (Turner et al. 2001, Strueder et al. 2001). Each of 
the MOS cameras provides an effective area of ^^600 cm^ at 
1.5 keV, and covers the energy range 0.2-10 keV; for the PN 
camera the corresponding values are, respectively, ^1400 cm^ 
and 0.15-15 keV. All the instruments were operated with the 
medium filter. The PN camera worked in Extended Full Frame 
mode in the first observations, with a CCD frame time of 200 
ms, and in standard Full Frame mode in the following ones. 
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Table 1. Summary of the observations 


Start time (UT) 
April 2000 

Off-axis angle 
arcmin 

Exposure 

ks 

MOS frame 

s 

PN frame 

ms 

7, 12h 35m 

0.2 

30 

2.6 

200 

9, 12h 22m 

10.3 

29 

2.6 

73 

11, 12h 26m 

10.4 

29 

2.6 

73 

15, 12h 26m 

10.2 

29 

2.6 

73 


with a frame time of 73 ms; both MOS cameras used the Full 
Frame mode in all the observations, with a frame time of 2.6 s. 

The first step of our data analysis was the event selection. 
For each of the four observations and for each of the three cam¬ 
eras we extracted all the events, with energy between 1 and 10 
keV, from a circular region with radius of 25" centered at the 
peak of the X-ray emission (RA = 18h 33m 33.8s, DEC = -10° 
34' 6 " (J2000)). This radius contains ~80% of the photons for 
a point source. We only considered events with pattern in the 
range 0-12 and 0-4 for, respectively, the MOS and the PA cam¬ 
era, in order to reject the non X-ray events due to cosmic rays 
and cosmetic defects. 

The times of the selected events were converted to the Solar 
System barycenter with the Reconstructed Orbit Files provided 
by the XMM Survey Science Center. These events were still 
tagged with discrete arrival times, corresponding to the readout 
times of the individual CCD frames: these arrival times were 
“randomized” by subtracting a random value between 0 and 
the relevant frame time from the original times. 

On these data we performed search for periodicities based 
on Fourier analysis (see, e.g., van der Klis 1989). A Fourier 
power spectrum was computed for each of the instruments and 
observations and examined for the presence of peaks above a 
threshold corresponding to a chance probability of 10 “^ of be¬ 
ing exceeded in the absence of a signal in a single spectrum. 
No significant peaks were found. 

To increase the detection sensitivity we repeated the same 
analysis on the sum of the individual spectra. To this aim, due 
to the three different CCD frame times (0.073, 0.2 and 2.6 s), 
we had to consider three different frequency ranges for the data 
analysis (Tab.2): 

- for ly < 0.1923 Hz = (2x2.6 s)“^, we summed the power 
spectra of all the observations and instruments. 

- for0.1923Hz< < 2.5 Hz = (2x0.2 s)“^, we summed the 

power spectra of only the PN data of the four observations. 

- for 2.5 Hz < ^ < 6.85 Hz = (2x0.073 s)“^, we summed 
the power spectra of the PN data of the three off-axis ob¬ 
servations. 

Again, no statistically significant peak was found in the 
summed power spectra. We also repeated the whole procedure 
for two distinct energy ranges, 1-3.5 keV and 3.5-lOkeV, again 
without any significant pulsed signal. 

To compute the upper limits on the source pulsed fraction, 
in the assumption of a sinusoidal pulse shape, we followed the 
procedure described in van der Klis (1989), taking into account 
the relevant correction pointed out by Vaughan et al. (1994). 
Finally, we had to correct the resulting upper limits to take into 


account the fraction of unpulsed flux due to the nebular emis¬ 
sion. We based this correction on the results of the Chandra 
observation (Slane et al. 2000), showing that at least 92% of 
the flux within 30" is of diffuse origin. Thus we obtained the 
upper limits (99.9 % confidence level) on the pulsed fraction 
of the neutron star emission given in Table 2, where we also 
report the number of counts used in the analysis. 

3. Discussion 

Our upper limits can be compared with those obtained by some 
recent works based on X-ray data of G21.5-0.9 . Note that 
the values reported by Warwick et al. (2001), 3.5% and 5.5% 
(respectively for MOSh-PN and PN only data), refer to the to¬ 
tal flux within an extraction radius of 8 ", including the nebular 
emission. Our results for the on-axis observation alone (the one 
used by these authors) are similar, but we reached a better sen¬ 
sitivity in the summed power spectra (corresponding to a factor 
~4 greater exposure time). 

For instrumental reasons, our search was limited to periods 
longer than 146 ms. Although it is clearly possible that the pul¬ 
sar in G21.5-0.9 has a shorter period, we note that Safi-Harb et 
al. ( 2001 ), from energetic considerations, estimated a period of 
P = 0.144 ( 745 /( 7 ^ 37 X 3 ))^/^ s (where I = 10"‘®745 g cm^ is the 
moment of inertia, E = lO^^i /37 erg s“^ is the spin-down en¬ 
ergy loss and r = T 3 x (3 kyr) is the pulsar characteristic age). 
For reasonable values of t and E, such a period falls in the 
range we could explore with EPIC. The same authors analyzed 
five data sets (total 75 ks) obtained with the Chandra HRC in¬ 
strument. They report an upper limit of 16%, without quoting 
the confidence level and, presumably, referring to the total flux 
within 2 ". 

The “canonical” picture of plerionic supernova remnants is 
based on young, energetic neutron stars with short rotation pe¬ 
riods, such as the Crab pulsar (P=33 ms) or the recently discov¬ 
ered pulsar in 3C 58 (P =66 ms, Murray et al. 2001). However, 
other results show that also sensitive searches for slower pul¬ 
sars are relevant: there are in fact relatively young pulsars with 
long periods. Besides the well known example of PSR B1509- 
58 in the SNR G320.4-01.2 (P=150 ms), other recent findings 
include the 325 ms pulsar in the SNR Kes 75 (Gotthelf et 
al. 2000), and PSR 111 19-6127 (P=407 ms, Pivovaroff et al. 
2001 ), which however does not have a bright synchrotron neb¬ 
ula. 
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Table 2. Upper limits on the pulsed fraction (99.9 % c.l.) 


Frequency range 

1-10 keV 

1-3.5 keV 

3.5-10 keV 

Hz 

% 

counts 

% 

counts 

% 

counts 

< 0.19 

6.2-7.5 

522142 

9.1-11.0 

300332 

12.4-13.5 

221810 

0.19 < 2.5 

15.1-17.1 

275543 

19.5-22.1 

154231 

24.2-27.3 

121312 

2.5 < 6.8 

24.6-30.1 

170434 

32-39.2 

98105 

33.3-40.8 

72329 
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